Capture into Rydberg states and momentum distributions of ionized electrons 
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The yield of neutral excited atoms and low-energy photoelectrons generated by the electron dy- 
namics in the combined Coulomb and laser field after tunneling is investigated. We present results of 
Monte-Carlo simulations built on the two-step semiclassical model, as well as analytic estimates and 
scaling relations for the population trapping into the Rydberg states. It is shown that mainly those 
electrons are captured into bound states of the neutral atom that due to their initial conditions (i) 
have moderate drift momentum imparted by the laser field and (ii) avoid strong interaction ( "hard" 
collision) with the ion. In addition, it is demonstrated that the channel of capture, when accounted 
for in semiclassical calculations, has a pronounced effect on the momentum distribution of electrons 
with small positive energy. For the parameters that we investigated its presence leads to a dip at 
zero momentum in the longitudinal momentum distribution of the ionized electrons. 

PACS numbers: 32.80Rm, 42.50Hz 



I. INTRODUCTION 

The interaction of strong laser radiation with atoms 
and molecules generates a variety of highly nonlinear 
phenomena. These include the extended photoelectron 
spectrum of above-threshold ionization (ATI) , the exces- 
sive yield of doubly and multiply charged ions, the gen- 
eration of very high harmonics of the driving field, and 
theyield of neutral atoms in excited Rydberg states (see 
[EH Sllll for reviews). The generally recognized trans- 
parent physical picture of the first three phenomena is 
based on the so-called direct ionization and rescattering 
scenario. According to it, electrons are promoted to the 
continuum via tunneling ionization and, while oscillating 
in the laser field, interact with their parent ions. Most 
of the ionized electrons, i.e. the "direct" electrons, are 
decelerated and deflected by the ion at small angles and 
finally contribute to the low-energy part of the photoelec- 
tron spectrum. In some rare cases, the ionized electron 
returns close to its parent ion and recombines emitting 
a high-frequency harmonic photon or kicks out one or 
several electrons, producing thereby a doubly or multi- 
ply charged ion. Yet another opportunity upon the close 
encounter is that the electron is elastically backscattered 
and thereafter acquires additional energy from the laser 
field. Such electrons contribute to the high-energy part of 
the photoelectron spectrum. The rescattering picture is 
explicitly realized in terms of complex "quantum orbits" 
or even simpler in terms of the two- and three-step semi- 
classical models. In the latter, real electron trajectories 
after ionization are calculated with Newton's equations. 

The interpretation of the afore-mentioned processes 
is largely classical, based on trajectories obtained from 
Newton's equations. Population trapping in Rydberg 
states resulting in an excited neutral atom at the end 
of the laser pulse is an exception because until recently 
its interpretation was purely quantum mechanical. The 
effect was considered as evidence of the existence of stabi- 
lization against ionization and studied theoretically and 



experimentally in many papers (see 0, [B[ and references 
therein). It is generally believed that there are two dif- 
ferent mechanisms of stabilization: adiabatic stabiliza- 
tion (stabilization in the Kramers— Henneberger regime) 
and interference stabilization. The adiabatic approach 
is appropriate when the field frequency lu exceeds the 
atomic ionization potential I so that the states in the 
discrete spectrum and the continuum are coupled by one- 
photon transitions. In experiments with optical lasers, 
such conditions were realized by initially populating Ry- 
dberg states. The physical mechanism behind interfer- 
ence stabilization is associated with multiple Raman-type 
transitions between the Rydberg levels and the common 
continuum. Destructive interference of the transition am- 
plitudes from these coherently populated states to the 
continuum suppresses ionization and, hence, serves to re- 
tain population in Rydberg states. 

The recent paper Q reports data on the yield of neu- 
tral excited He* atoms and singly charged ions He + from 
a gaseous target irradiated by 30-fs Titanium-Sapphire 
(Ti:Sa) laser pulses with intensities up to 10 15 W/cm 2 . 
The data, corrected for the radiative decay of neutral 
excited atoms on the way to the detector, fairly well 
agree with the results of extensive quantum-mechanical 
calculations. The latter were accomplished in two ways: 
by approximately solving the two-electron TDSE for he- 
lium exposed to a laser field and in single-active-electron 
(SAE) approximation. The excited neutral atom yield is 
roughly 20% of the ion yield at lower intensities decreas- 
ing to about 10% at higher intensities. 

In addition, a two-step semiclassical model was em- 
ployed for the first time for providing a simple picture 
of the dynamics resulting in the population of Rydberg 
states. Monte-Carlo calculations were performed with 
classical trajectories in the combined Coulomb potential 
and electric laser field. Those electrons are captured into 
bound states of the neutral atom that have negative to- 
tal (kinetic plus Coulomb) energy at the end of the laser 
pulse. The energy distribution in the subset of bound 
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trajectories was converted into the distribution in an ef- 
fective principal quantum number (E = — 1 /2n 2 ) and the 
latter was found to be in qualitative and even quanti- 
tative agreement with the one deduced from quantum 
calculations. 

The implication of this study is that in Monte-Carlo 
simulations with classical trajectories in the laser field 
and the ionic Coulomb field not all electrons that tun- 
neled out at some time are actually ionized after the end 
of the laser pulse. Some may be captured back into bound 
states. This outcome of the electron motion in the con- 
tinuum was dubbed frustrated tunneling ionization (FTI) 
in Transient trapping in Rydberg states during the 
laser pulse was mentioned earlier in [7[. However, to the 
best of our knowledge, there are no comments on the for- 
mation of excited neutral atoms at the end of the pulse in 
the papers where Monte-Carlo simulation of classical tra- 
jectories after tunneling was used to model the effect of 
Coulomb focusing on the rescattering processes [7|, [8|, l2( 
or to investigate the low-energy photoelectron spectra 

In this paper, we further exploit the semiclassical 
model: (i) to get deeper insight into the mechanism re- 
sponsible for the yield of neutral exited atoms; (ii) to find 
out the scaling of this yield with the laser-atom parame- 
ters; and (iii) to understand the impact of the population 
trapping on the momentum distributions of the ionized 
electrons. 



II. MODEL 

In this section we will sketch our simulation technique 
with special emphasis on the details that are essential for 
the following. 

Let us consider a neutral atom with the ionization po- 
tential /, irradiated by a short laser pulse with the du- 
ration tl = (2Tr/to)n p and sine-square envelope, linearly 
polarized along the x axis: 

F(t) = F sin 2 (t/r L ) cos (ut) , (1) 

where n p is the number of cycles and u> is the carrier 
frequency. At time to the bound electron tunnels out of 
the atom quasistatically with the rate 

l \ 3F(t a ) [ 2/J J w 

Here I and F a = (2/) 3 / 2 are the ionization potential 
and the characteristic atomic force of the bound state 
(I = 1/2, F a = 1 for the ground state of hydrogen). 
The rate © is obtained from the standard strong-field 
ionization amplitude [H, [l4| calculated by the saddle- 
point method in the tunneling regime, when the Keldysh 



parameter 7 = uj\/2J/F <C 1. Details of the calculations 
are presented elsewhere (see [HI EH)- Atomic units e = 
m = h = 1 are used. 

At the tunnel exit xq = I/F(to) the electron has 
zero longitudinal and nonzero initial transverse velocity 
vq. The released electron moves along a classical trajec- 
tory under the action of the laser field and the Coulomb 
force of its residual ion. The initial conditions unam- 
biguously determine its position r (t) = r(to,vo,t) and 
velocity q(t) = q(to,i>0)i) at any time instant t. In the 
limit of large i, one obtains the asymptotic momentum 
p = q (to, vo,t — > 00), which is identical with the momen- 
tum measured by the detector. During the whole laser 
pulse < to < tl, an ensemble of trajectories is launched 
corresponding to different pairs (to, i>o) weighted with the 
probability ([2]). No interference effects are accounted for. 
In our simulations we employed about 1.5 • 10 6 trajecto- 
ries regardless of the intensity. By solving the respective 
Newton equations one obtains r (t) and q (t) for each tra- 
jectory. The ratio of the number of trajectories in a par- 
ticular bin of momentum space to the total number of 
trajectories is the momentum distribution at given time 
t normalized to unity. It evolves in time and in the limit 
t — ► 00 can be compared with the experimental data. 

The picture described in the previous paragraph is not 
complete, however. Indeed, if the Coulomb field of the 
nucleus is included into the Newton equations then, for 
some initial conditions (to, vq), the electron's momentum 
simply does not converge to any constant value but keeps 
oscillating in the limit t — > 00. Such trajectories corre- 
spond to those electrons whose total energy after the end 
of the laser pulse is negative 

E(r L ) = ^---<0, (3) 
2 r 

so that they never reach a detector. Following Q we 
interpret these electrons as trapped in Rydberg states 
(captured electrons). Here q = q_(v , t , tl) and r = 
t(i>o, to,TL,) are the electron velocity and its distance from 
the ion at the end of the laser pulse, respectively. Starting 
with the time tl the electron moves only in the Coulomb 
field, and its energy as well as the angular momentum 
are conserved. Thus, since this time the whole ensemble 
of trajectories can be subdivided into two subsets. The 
subset of trajectories satisfying condition ([3]) specifies the 
distribution of the bound electrons over the principal and 
the angular quantum numbers n and I. The subset of tra- 
jectories with positive total energy determines the distri- 
bution of ionized electrons. 



III. RESULTS AND DISCUSSION 

The results of our simulations are presented in Figs. [1] 
H and [31 

The energy distribution of the electrons at the time tl 
is shown in Fig. [T] It can be seen from the figure that 
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FIG. 1: The distribution of the total energy E = q 2 /2 - 1/r of 
electrons tunneling out of an H-atom (/ = 13.6eV) at the end 
of a laser pulse with a duration of n v — 5 cycles, frequency 
u = 0.05 a.u. and an intensity of 3.5 • 10 14 W/cm 2 . The 
number of electrons with negative energy (shown by green 
bars) is estimated at about 15%, which is consistent with [f|. 



the number of captured electrons corresponds to approx- 
imately 15% of the number of ions, which agrees with the 
results obtained in The distribution of the trapped 
electrons over the principal quantum number is depicted 
in Fig. O This distribution is similar to those evaluated 
in 0] by the solution of the time-dependent Schrodinger 
equation and by Monte-Carlo simulations: the maximum 
of the distribution corresponds to n = 6. 

Finally, Fig. [3] exhibits the areas of the plane (to,Vo) 
of initial conditions that lead to capture of the electron 
into Rydberg states. The profile of the laser pulse is also 
shown in the figure. Figure [3] reveals two main features. 
First, the capture does not take place for any initial veloc- 
ity. It occurs only when \vq\ is less than some maximum 
value (approximately 0.3 for the parameters of Fig. [3]). 
But this is not yet a sufficient condition. Such a con- 
dition will be formulated in the next section, in which 
Fig. [3J is explained and all necessary estimates are done. 



A. Mechanism and scalings of population trapping 

To understand the results of Fig. [3J one should keep in 
mind that the momentum (velocity) q and the distance 
r in the capture condition ([3]) are actually correlated. 
The larger the momentum, the larger is the electron dis- 
tance and the smaller is the potential energy and, hence, 
the electron cannot be trapped. With the equality sign, 
Eq. ([3J sets an upper limit for q. Regarding smaller val- 
ues for q, the electron can have a small drift momentum 
at the end of the laser pulse only under certain condi- 



tions. The trajectories of electrons ionized near the field 
extrema with initially small drift momentum along the 
laser field and with small transverse velocity are strongly 
affected by the attraction of the Coulomb field. While 
such an electron oscillates in the laser field, a single 
"hard" collision or a sequence of "soft" collisions Q are 
likely to increase its drift momentum so that, with very 
high probability, it is not a candidate for capture any- 
more. On closer inspection, this holds for electrons that 
start their orbit after an extremum of the field. In this 
case, the drift momentum imparted by the laser field is 
directed towards the ion. The net effect of the laser field 
is that these electrons are driven back to the ion at least 
once if not many times. An example of such a trajectory 
is shown in Fig. H]by the solid line. The Coulomb force 
still further increases the drift momentum in the direc- 
tion towards the ion. Inevitably electrons with large drift 
momentum return back to the ion and scatter depend- 
ing on the magnitude of their impact parameter. The 
opportunity to avoid strong interaction with the ion and 
have sufficiently small drift momentum at the end of the 
laser pulse exists for electrons which start in the contin- 
uum before the field extremum. Without the Coulomb 
field, while oscillating in the laser field they drift away 
from their parent ion and never return to its vicinity 
(the trajectory shown by the dashed line in Fig. 2]). The 
Coulomb attraction reduces the drift momentum without 
reversing its direction, so that the trajectory remains "on 
the average outgoing" . This classification of trajectories 
is illustrated by Fig. 0J 

To support this general conclusion, let us consider the 
contribution to the capture process from a single half cy- 
cle of the laser field, for example, from the central maxi- 
mum of the field ([!]) . With time counted from this max- 
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FIG. 2: The distribution of the captured electrons over the 
principal quantum number n. The parameters are those of 
Fig. 1. 
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FIG. 3: (a) The 2D distribution of the captured electrons 
over the time of tunneling to and initial velocity vo for the 
parameters of Fig. 1. The laser field {JTJ is shown by a blue 
line, (b) The same as in (a), but for the dominant maximum 
of the pulse. Dark means high electron yield. 



imum and for phases within the interval — n/2 < uot < 
7r/2, the field and its vector potential are simplified as 
F(t) = Fcos(wt) and A (t) = - (F/u) sin {uxt). The 
momentum at the end of the laser pulse of an electron 
tunneling at time to with nonzero transverse velocity can 
be written as 

q= q(t ,v ,T L ) = A (t ) + v + p c (t ,v ). (4) 

At t > tl the momentum (Q} evolves in the Coulomb 
field. The trajectory starts at the tunnel exit xq — 
I/F (to) and the field F (to) accelerates the electron away 
from it. For the laser-atom parameters of interest, the 
Coulomb force is small in comparison to the laser field 
already at the tunnel exit and rapidly decreases when 
the electron departs from the ion. For example, with the 
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FIG. 4: Two electron trajectories in the field F(t) = 
Fo cos(wf) corresponding to two instants of ionization before 
and after the field maximum: unto = —0.1 (dashed curve) and 
tvto = 0.1 (solid curve) for the parameters of Fig[T] The coor- 
dinate is scaled by the amplitude Ao = Fo/lo 2 of the electron 
oscillation in the laser field. 



parameters of Figs. 1 and [4] and not too far from the 
field maximum one has xo ~ 5 a.u., while the oscillation 
amplitude A = F /oj 2 is about 40 a.u. This allows us to 
evaluate the contribution of the Coulomb force in Eq. ^ 



PC (to, v ) 



to ' L 



r\ (t) 



(5) 



as the integral along the initial part of the electron tra- 
jectory 

r L (t) = jz + ±F (to) (t - t ) 2 ,v (t- to) | , 

and to extend the upper integration limit to infinity [l5| . 
Evaluation of the integral ([5]) for not very large transverse 
velocities, vo < V2I, gives the drift momentum (j4]) in the 
plane defined by the directions of the field (x axis) and 
of the transverse initial velocity (_L axis): 



F . Fcos(wto) 
q x = --sin Mo)-* (2j)3/2 , 

9 \F(t a )\ 



(G) 
(7) 



Except for the numerical factor, the Coulomb contribu- 
tion to the drift momentum © is just the product of the 
Coulomb force at the tunnel exit, — l/xfj, and the time 
interval At = yj2xo/F (to) over which the electron-ion 
distance doubles. Not surprisingly, this time interval co- 
incides with the time of flight under the potential barrier, 
Xo/v2I, introduced by Keldysh. 

For the next half period, the laser field changes its 
direction and pulls the outgoing electron back to the ion. 
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But, for q x > 0, the minimum distance in the x direction 
between the electron and the ion will be larger than the 
tunnel exit xq by the drift displacement q x /uj. So, if 
this displacement is of the order of or larger than Xq, 
the Coulomb force upon the closest approach is less than 
the one at the time of tunneling, and its effect can be 
neglected. In such cases, the drift momentum given by 
© and {7} remains practically constant during the rest of 
the laser pulse. This allows us to estimate the electron- 
ion distance at the end of the laser pulse in Eq. ([3]) as 



= <7 (t l - t ) . 



(8) 



where q = v^x + 9± ■ With the momentum and the dis- 
tance at the end of the laser pulse known as functions of 
to and i>o, the condition ([3]) determines the subspace of 
the initial parameters that result in capture to the bound 
states: 



q < 



T~L — to 



1/3 



(9) 



It is immediately obvious that the number of trapped 
electrons decreases with increasing time that the electron 
spends in the laser field after ionization. This explains 
the difference in the number of electrons captured from 
identical laser half-periods before and after the field max- 
imum, which is seen in Fig. [3Ji. 

Next, we try to obtain a rough analytical estimate for 
the relative yield of neutral excited atoms with respect 
to singly charged ions, i.e. the ratio N* /N + , in a short 
laser pulse. The numbers N* and N + can be evaluated 
by integrating the rate ([2]) over the respective parts of 
the (to, ^o)-plane shown in Fig. [3K To a fair approxi- 
mation, the integration can be performed over the initial 
conditions related to the central maximum of the field. 
Furthermore, applying the theorem of the mean to both 
integrals, we assume that the rates ([2]) do not differ essen- 
tially and cancel out in their ratio, which thereby is esti- 
mated as N*/N+ m E*/E+. Here E* and E+ are the ar- 
eas of the (to, vq) plane effective for population trapping 
and real ionization, respectively. For ions from the cen- 
tral maximum of the field we have E + w n x \p2J^I F/F a , 
where the second factor is the width of the rate in the 
initial velocity. 

For the neutrals, E* is just the area of the crescent- 
shaped region in Fig. [3}d. First, we find the interval of 
start times that are favorable for capture for vo = 0. Its 
left end io< is found by equating the two sides of Eq. 
([9]) under the assumption that we deal with an "on the 
average outgoing trajectory" , i.e. q x > 0. The resulting 
equation is 



sin (ojto) — A cos (coto) 



F 



TL - to 



1/3 



(10) 



where A = nu;/ (2I) 3 ^ 2 . In the middle of the laser pulse, 
tl — to ~ Tl/2, and, typically, the r.h.s. of Eq. (flU)) is 



small in comparison with unity. To first order in coto < 1 , 
the solution within the interval — n/2 < ujto < 7r/2 is 



uto = ooto< 



F \t l 



1/3 



(11) 



The momentum © decreases when the start time to > 
to< moves away from the left end of the favorable in- 
terval. If it is recalled that trajectories with very small 
momenta are strongly perturbed by the Coulomb field 
and do not contribute to the process of capture, then 
the upper end of the window can be estimated from the 
equation q x (to>,vo = 0,T£, = 0). For A < 1, one easily 
finds 



ujt 



o> 



-A. 



(12) 



For the parameters of Fig. [31 the predictions of Eqs. (fTTj) 
and (TT2"1) . i.e. wto< ~ —15° and ujto> ~ —9°, respectively, 
agree nicely with the results of the Monte-Carlo simula- 
tions presented in the lower panel of this figure. 

It can be seen that the transverse width of the crescent- 
like shape in Fig. Ob, i.e. the maximum transverse veloc- 
ity allowing for the capture of electrons into the bound 
states, is more or less the same for start times within the 
favorable window. To estimate this width, we consider 
the condition ([9]) for the case when q x (to>, Vq, t£) = 
and for times in the middle of the laser pulse so that 
tl — to ~ tx/2. With account of 0, it immediately 
takes the form \v$\ < v maK where 



1/3 



/ 1-2 



F 



(2iy 



(13) 



Actually, in (fl3|) we have neglected terms of the order 
of A 2 -c 1. In Fig. O the transverse velocity does not 
exceed 0.30 a.u. whereas, for the same parameters, one 
has from Eq. (fT5j) w max = 0.23 a.u. Again, the agreement 
is quite reasonable. 

By roughly estimating the capture area in Fig. [3] as 
E* w oj (to> — £()<) x 2t> max and dividing it by the area 
E + effective for ionization (see above) we have for the 
relative yield of neutral excited atoms with respect to 
singly charged ions 



N* 
N+ 



F 3/2 2/3 



1 - 2- 



F 



(2iy 



(14) 



The scaling with the laser parameters given by Eq. (|14p 
was correlated with the results of the Monte-Carlo simu- 
lations. According to the latter, when the field strength 
varies from F = 0.1 a.u. to F = 0.07 a.u. the ratio 
N*/N+ rises from 0.18 to 0.41, i.e. by a factor of 2.3, 
while Eq. (fT4"]l predicts a factor of 1.7. Agreement is even 
better with respect to pulse duration. When the pulse 
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FIG. 5: The longitudinal momentum distribution evaluated 
by Monte-Carlo simulations at the end of the laser pulse at 
time tl- Dotted blue curve: All launched trajectories, includ- 
ing those with E(jl) < 0, are accounted for. Solid red curve: 
The trajectories with negative energies are excluded. Dashed 
green curve: The result of simulations with the Coulomb field 
entirely turned off. The parameters are those of Fig. 1. 



duration is doubled from 5 to 10 cycles, then for the field 
strength of F — 0.1 a.u. the relative yield drops from 0.18 
to 0.11, i.e. by a factor of 1.64, which is practically equal 
to the factor of 2 2 / 3 expected from (fl"4")l . The general 
tendency of an increasing percentage of bound electrons 
for shorter pulses and lower intensities was observed in 
the Monte-Carlo simulations of Ref. • 

Experimentally, a strong decrease of the yield of neu- 
tral excited atoms with increasing pulse duration was re- 
ported in an early paper [ljj. When the scaling of Eq. 
(|14|) is applied to analyze ionization of He atoms with 
their high ionization potential, the very last factor can 
be dropped. Then it follows from (fT4]) that an-order-of- 
magnitude increase of the laser intensity will reduce the 
relative yield by a factor of 5 to 6. The relative yield 
reported in Ref. Q decreases over such a range by ap- 
proximately a factor of 2. The origin of this discrepancy 
is not presently clear. 

According to [6], the maximum of the n-distribution 
scales with y/F /u. The derivation of this dependence 
presented there was based on the heuristic assumption 
that the radial expectation value of the final Rydberg 
atom, (r) oc n 2 , is proportional to the amplitude of elec- 
tron oscillation in the laser field, Aq — F/uo 2 . Our idea 
is that the quantum-mechanical mean distance of a cap- 
tured orbit is determined by the electron distance at 
the end of the laser pulse, (r) w r (to). From Monte-Carlo 
simulations with the parameters of Fig. [TJ we found that 
the distribution of the distances of trapped electrons (not 
shown here) has a maximum around r max = 60 a.u. The 




FIG. 6: The longitudinal momentum distributions, calcu- 
lated taking into consideration all launched trajectories, at 
t = 2tl (dotted blue curve), t — 3tl (dashed green curve), 
t = 4tl (dash-dotted magenta curve). The red curve shows 
the asymptotic distribution, which was evaluated by exclud- 
ing the trajectories with E(tl) < 0, as discussed in the text. 

relation r max = (r) — 3n 2 /2 then gives n max ss 6, which 
agrees nicely with the results of Fig. [2l The distance ([8]) 
indeed becomes proportional to the oscillation amplitude 
if one factors out the dimensional factor F/lo from the 
momentum and converts time into the field phase. 

B. The capture into Rydberg states and the 
distribution of the ionized electrons 

The fact that a substantial part of the tunneled elec- 
trons end up in bound states with E < after the end 
of the pulse should affect the energy- angular spectrum of 
direct ionization, particularly its low-energy part. This 
raises the question of a possible relation between the cap- 
ture into Rydberg states and the dip (minimum) at p = 
in the momentum distribution of ionized electrons along 
the laser polarization, which was for the first time ob- 
served in [l8j. The origin of this dip has been the object 
of intense investigations [H [Q, EI US HI, SU . In this 
section we analyze momentum distributions of the ion- 
ized electrons, calculated taking into account the capture 
into bound states. 

First, however, we address the photoelectron momen- 
tum distribution at the end of the pulse (not at infinity) 
(see Fig. [5]), calculated with the assumption that all elec- 
trons are ionized, regardless of whether the total energy 
E(tl) [Eq. ([3])] is positive or negative [Io| . Note, that in 
the presence of the Coulomb field such a distribution is 
not an observable, since the momentum is not conserved. 
However, the momentum distribution at the end of the 
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FIG. 7: (a) The 2D momentum distribution for the parame- 
ters of Fig. 1 resulting from Monte-Carlo simulations at the 
end of the laser pulse, where the trajectories with negative 
total energy E = q 2 /2 - 1/r < are discarded, (b) The 
same distribution in the asymptotic limit after its subsequent 
evolution in the Coulomb field. 



pulse is instructive for understanding the influence of the 
Coulomb force on the electron dynamics. The distribu- 
tion evaluated without taking the Coulomb field of the 
atomic residual into account is also shown in the figure. 
Figure [5] demonstrates that the Coulomb force severely 
modifies the central part of the longitudinal momentum 
distribution. Incidentally, the asymmetry of the dip is 
due to the small (5 optical cycles) duration of the laser 
pulse. 

Figure [6] illustrates the evolution of the distribution 
of p x calculated considering all launched trajectories, in- 
cluding those with E(tl) < at the end of the pulse, 
under the action of the Coulomb force. One observes 
that the dip in the longitudinal distribution is gradually 
filled in and a narrow maximum develops in its place. 
Our simulations show that the dip in this approach also 
vanishes with increasing laser intensity and pulse dura- 
tion. We should focus our attention on the fact that 
the latter result coincides with a conclusion of Ref . [ll| , 
which, however, concerned the asymptotic distribution of 
the ionized electrons. It is very important to note that 
the asymptotic distribution, by which we understand the 
distribution to be measured by a detector far away from 
the ion, has to be calculated by excluding the trajecto- 
ries with E{t l ) < 0. Since the latter are trapped into 
Rydberg states they will never reach the detector. If nev- 
ertheless they are included into the asymptotic distribu- 



tion, their momenta will not converge to definite values as 
time goes to infinity or, depending on the numerical ac- 
curacy, they will contribute to a spurious accumulation 
of events with very low momenta. Indeed, the asymp- 
totic distribution obtained by excluding the trajectories 
with negative energy (the solid red curve in Fig. [SJ has a 
distinct minimum at zero longitudinal momentum. 

Let us summarize the results shown in Figs. [5] and [5] 
A dip in the longitudinal momentum distribution gener- 
ated by a relatively short laser pulse is already seen at 
the end of the pulse. A comparison of the two curves 
calculated with and without the Coulomb field suggests 
that the Coulomb field is responsible for the formation of 
this dip [Hi EI- However, Fig. [5] shows that this dip is 
washed out if we allow the distribution to continue devel- 
oping under the action of the Coulomb field after the end 
of the laser field. Our definition of the asymptotic distri- 
bution discards the contribution of those electrons with 
negative energy at the end of the pulse. Remarkably, the 
asymptotic distribution again exhibits a dip. The obvi- 
ous conclusion is that the mechanism responsible for the 
formation of the dip at zero momentum is closely related 
to the one that governs the electron capture into Rydberg 
states. Both effects are caused by the Coulomb field and 
occur in adjacent energy intervals. We notice in passing 
that the observed structure also survives focal averaging, 
which is inevitable in real experimental conditions. 

Hence, the calculation of the momentum distribution 
of ionized electrons must take into account the possible 
population of Rydberg orbits with negative total energy. 
The two-dimensional momentum distribution, which dis- 
regards electrons with E{t) < at t — tl, is depicted in 
Fig. [7(a). One can see from the figure that the central 
part of the distribution is completely unpopulated. Sim- 
ulations for t > tl show that the subsequent motion in 
the Coulomb field partly fills in this hole, see below. 

The asymptotic distribution at t — > oo can be obtained 
without calculations up to a large time instant i on d, which 
guarantees the asymptotic regime for the vast majority 
of launched trajectories. Indeed, the electron momen- 
tum q(to)Wo,Ti) and its position vector r(to, Vq, tl) at 
the end of the laser pulse uniquely determine the subse- 
quent trajectory in the Coulomb field of the parent ion. 
Standard formulas of classical mechanics for hyperbolic 
motion (see, for example (23j). allow straightforward an- 
alytical evaluation of the electronic asymptotic momen- 
tum 

p(A x M) - A 



V 



(15) 



1 + p 2 M 2 ' 

where M = r x q and A = qxM-r/r are the conserved 
angular momentum and Runge-Lenz vector, respectively. 
The absolute value p of the asymptotic momentum can 
be determined from the energy conservation law: 



(16) 



q 2 1 p 2 
T ~ r ~ T' 

It should be noticed that Eq. (|16p can also be used to 
explain the origin of the unpopulated area in Fig. [7Ja): 
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FIG. 8: Asymptotic momentum distributions along laser po- 
larization for different laser intensities (in W/cm 2 ) for the 
parameters of Fig.l. 

In fact, at the end of the laser pulse electrons can be 
found at various finite distances from the ion. The pos- 
itivity of the total energy imposes a lower limit on the 
value of the electron momentum at that instant of time. 
In contrast, the electron does not feel the Coulomb field 
at t — » oo, and its momentum can be arbitrarily small. 

The distribution evaluated by this approach is shown 
in Fig. [7](b) . Together with Fig. [8] it demonstrates 
that the dip is not smeared out by the evolution in the 
Coulomb field after the end of the laser pulse and, more- 
over, its shape is not very sensitive to the laser intensity. 
The same holds for the behavior of the distribution with 
increasing pulse duration. 

Thus, depletion of the low-energy part of the ioniza- 
tion spectrum due to the capture into Rydberg states is 
a possible mechanism for the formation of a dip in the 
momentum distribution along the polarization direction. 

IV. CONCLUSIONS 

Based on the two-step semiclassical model, we have in- 
vestigated the electron motion after tunneling in a strong 



laser field considering the interaction with the Coulomb 
field of the parent ion. The focus of attention is on the 
mechanism governing the yield of neutral excited atoms 
from a gaseous target irradiated by a short laser pulse. A 
simple physical picture emerges from Monte-Carlo simu- 
lations with classical trajectories after tunneling. With 
account of the Coulomb field, not all tunneled electrons 
are actually ionized since some of them are captured back 
into bound atomic states after the end of the laser pulse. 
Obviously, an electron remains bound if its total energy 
at the end of the laser pulse is negative. Furthermore, 
we show that for an electron to be captured it (i) has to 
have moderate drift momentum at the time of tunneling 
and (ii) has to avoid strong interaction ( "hard" collision) 
with the ion. Such trajectories do exist when electrons 
are released within a narrow time window before the ex- 
trema of the oscillating laser field having a not too large 
transverse velocity Vo- Understanding these features has 
allowed us to derive the scaling of the ratio of neutral 
excited atoms and singly charged ions with the laser- 
atom parameters in analytic form. Its prediction, that 
the percentage of neutrals increases with decreasing in- 
tensity or pulse duration, fits into the general physical 
picture: In both cases, the electron drifts away from the 
ion by a smaller distance during the laser pulse and ex- 
periences a stronger Coulomb attraction. This distance 
is smaller at lower intensity because of a reduced drift 
momentum imparted by the laser field and, in the other 
case of a shorter pulse, because of a smaller travel time. 
In addition, it has been demonstrated that the channel of 
capture, when accounted for in semiclassical calculations, 
has a pronounced effect on the momentum distribution 
of electrons with small positive energy. In particular, it is 
correlated with the formation of a "dip" at zero momen- 
tum in the distribution of the longitudinal momentum. 
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